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Three ynelhods of approximation are described and are used to separate the primary 
twilight hrighiness from the observed brightness. Photoelectric observaMons obtained at 
the Majaky Observatory in 1968-69 {Zaginayilo) are combined with observaMons from 
a balloon and from Pic-du-Midi Observatory {Link et al.j 1967) to derive the atmos- 
pheric scattering phase functioiis at 0.37 jx and 0.58 p as a function of height. 

Comparison of these data with data for a Rayleigh atmosphere provide information 
on the optical properties of dust in the upper atmosphere. 


T he observed brightness Bobs of the twilight 
sky may be considered as the sum of the 
primary brightness Bi and the secondary one 

Bq\)b = B\-\-B2 ( 1 ) 

For investigation of the optical properties 
of the Earth^s atmosphere it is necessary to 
separate the primary brightness Bi from the 
observed brightness Bobs- It is a difficult problem. 
As there is no widely adopted method now for 
separating the primary brightness, we have 
determined this brightness by using three ap- 
proximation methods. 

The first method was propost^d by V. G. 
Fessenkov (1968a) and utilizes observations in 
two points of the solar meridian with zenith 
distances By this method the brightness 

of the secondary twilight in the observation 
point is determined from observe ations in the other 
point (th(^ base point) lying in the Earth's 
shadow and having the zenith distance ^ = 70"". 

This method may be extended to points of the 
sky with zenith distances different from ^ = 70°, 
Values of coefficient k arc listed in table 1. The 


coefficient k is used for computations of the 
secondary brightness B^iz) in the observation 
point from the secondary brightness 52^(70°) 
in the base point z = 70°, A = 180° by the formula 
B2(2)=A:B2'(70°). 

In the second method, proposed by Divari 
(1971), we use another coefficient ki which is 
defined by the ratio Bi{Bi + B 2 ) ^ki. This 
coefficient ki may also be found from the com- 
putations for the standard atmosphere. Values 
of this coefficient arc listed in table 2. 

These values found according to the data of 
Divari and Plotnikova (1965), for corresponding 
wavelengtlis and atmospheric transmission co- 
efficients may be considered only as a first 
approximation. 

The third method for the determination of 
the primary brightness is V. G. Fessenkov^s 
method of obsiTvations of the point of the sky 
mth zenith distance z = 80° and azimuth A = 30° 
(from the Sun azimuth). This method requires 
measurement of intensity, d('gree and angle of 
polarization of the light of the twilight sky in 
the quoted point (Fessenkov, 1968b). 
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Table 1. — Values of Coefficient k =Bi{z) /Bf {70°) used for A=0° at Wavelengths \=0.37 and \ = 0.5<9 m With 
Corresponding Atmospheric Transjnission Coefficients p=0.59 and p=0.89 {Base Point z=70°, \ = 180°) 



X = 0 . 37 m , p = 0.59 

X = 0-58 M, p = 0.89 
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0.684 
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Table 2. — Values of Coefficient Ki =Bi/(Bi-}-B 2 ) Use 1 for the Wavelengths \=0.37 and \=0.58 p {Transmission Coefficients 
p=0.S9 and p —0,89) and Computed According to the Data of Divari and Plotnikova {1965) for A=0° 


0 \ 

X = 0.37 p = 0.59 

X = 0.58 li , p = 0,89 

4 ° 

6 ° 

10 ° 

14 ° 

4 ° 

6 ° 

10 ° 

14 ° 
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0.66 

0.43 
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Having the primary briglitness Bi as tho func- 
tion of the Sun depression g{g^z—00°, where z 
is the zenith distance of the Sun) or of the (effec- 
tive height of the twilight ray it is possibles to 
determine the scattc^ring coefficient a{0, h) of 
tlie atmosphere as a function of the scattering 
angle 6 and th(' hc'ight h above the surface of the 
Eart^h. This procedure is described by Divari 
(1971). 

Knowing a(6y A), we can find tho volume 
scattering coefficient 2(A) (or extinction coeffi- 


cient) according to 

i:{h)=Ja{e,h)dw ( 2 ) 

where dw is an element of solid angle and inte- 
gration is over all the angles. 

For the purpose of comparing the three fore- 
going methods of separating primary brightness 
Bi and determining the scattering coefficients, 
we have used the photoelectric observations of 
twilight sky made by Y. I. Zaginayilo in 1968- 
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1969 at the observatory Majaky near Odessa 
(USSR). 

Figure 1 shows the observed brightness Bobs in 
the solar meridian at the points z=70° and 
2 =—60° (the minus sign indicates that the 
observed point is on the antisolar side of the Sun 
meridian) and the primary brightnesses Bi'‘> 
and Bi‘*> found by the two methods above. As 
may be seen, the Bi<') and Bi<“) at the point 
2 = 70° are nearly equal in the interval 2° to 8° of 
the Sun depression. For g>8°, the second method 
leads to lower values than the first one (Bi^^' < 
Bi'”). At the point z=-60° the deviation of 
Bi«” from Bi'*’ begins earlier, approximately at 
(7=4°. 

Figure 2 shows the primary brightnesses 
Bi<*^ and Bi<“^ at the wavelength X = 0.58;i ob- 
tained from the observed brightness Bobs at 
the point z = 30° of the Sun meridian. In the 



Pigtjuk 1. — T.fOg of the observed lirightiiess I^ohs Jit the 
points Z = 70® (crosses) and Z--60'' of the solar 
meridian and primary briglitnesses (open circles) 
and (black dots) for the same point of the sky, 
plotted against Qq, the solar depression angle. 


same figure we have shown the brightness Bi(25) 
(X = 0.53/i) which was obtained by F. Link, 
L. Ncuzil, and I. Zacharov (1967) by means of 
balloon observations at the height of 25 km 
above the Earth's surface. By crosses we have 
shown the ob serv ed brightnesses, obtained by 
Link et al. (1967) from the observation (at 
the wavelength X = 0.53 m) at the Pic-du-Midi 
observatory (.ff=2.9 km). It seems from com- 
paring the data presented in figure 2 that 
gives a more realistic approximation of primary 
brightnesses than 

In figure 3 we have shotvn the phase functions 
a (6) for X0.37 /u and X0.58 m obtained from the 
a(df h) by averaging in the height interval 30 to 
130 km. The vertical bars indicate the standard 
errors. The same figure shows the Rayleigh 
scattering function ^l + cos^0. As seen from 
figure 3, the intensity of scattered light increases 
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Figure 2. — Log of the observed brightness Bob# and pri- 
mary brightnesscis Bi^) and Bi^®^ at wavelengths X 
0.58 M from observations at the point Z = 30*’ of the 
sun meridian, Bi(25) is primary brightness according 
to Link et al. (1967) from balloon observation. Crosses 
are their ohseryations.of brightness at the Pic-du-Midi. 
Abscissae are solar depression angles qq. 





284 


EVOLUTIONARY AND PHYSICAL PROPERTIES OF METEOROIDS 


^vith the decrease of angle 6 steeper than for 
Rayleigh scattering. It seems that it is the mani- 
festation of the influence of aerosols on the 
optical properties of atmosphere. It may be 
pointed out that the increase of scattering 
intensity wdth the decrease of the scattering 
angle 0 is flatter at X0.37 m than at X0.58 m- 

Figure 3 also shows th(^ aerosol phase function 
(Taid) obtained from the atmospheric scattering 
coefficient a{6) by subtracting the Rayleigh 
scattering coefficient calculated for stand- 

ard atmosphere CIRA, 1965. In the same figure 
we show the aerosol phase function obtained by 
E. V. Pyaskovskaya-Fessenkova (1970) from 
observations of the day skylight. As may be 
seen, there is a good agreement between the 
phase functions obtained from twilight and 
da>dime observ^ations. 

The scattering coefficient <r(d) was defined 
from observations in the sun mc^ridian for the 
interval of scattering angles 0=25° to 155°. Since 
for the determination of the extinction coefficient 
S it is necessary to integrate (x{6) over the in- 
terval $ from ^ — 0° to 0—180°, we extrapolate 
this function from 0=25° to 0 = 0° and from 



Figuke 3— Atmospheric phase function < 737 ( 9 ) for X 0.37 n 
and <73g(9) for X 0.58 n and aerosol phase function 
<ra{$) for X 0.37 fi. P-F is aerosol phase function ac- 
cording to Pyaskovskaya-Fessenkova (1970) for X 
0.546 Ai. <7fl(9) is Rayleigh phase function for X 0.37 fi. 


0=155° to 0=180°. The extrapolation from 
0=25° to 0 = 0° was made by taking into account 
the results found by G. Newkirk and J. A. Eddy 
(1963) for the small angles from the stratospheric 
coronagraph observations. The extrapolation in 
the angle interval 160° to 180° is made with 
sufficient certainty and the error of this extra- 
polation according to our estimations can reach 
a few percent of the value of S. 

Figure 4 shows the extinction coefficients 
2 ( 1 )^ 2 ^ 2 ) 2J<3) obtained for the primary 



FiGimE 4. — Extinction coefficient 2^)^ and 2^®^ ob- 
tained from the primary brightnesses and 

at the wavelength X 0.37 fx. Point X is the value 
according to Mikirov (1963). Dotted line portrays 
extinction coefficient according to Fiocco and Colombo 
(1964), 
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Figure 5. — Ratios 2/2ij. Curve 1 is ratio obtained from 
in the Sun meridian. Curve 2 is ratio obtained 
from in the Sun meridian. Cur\-^c 3 is ratio ob- 
tained from at the point Z =80°, A =30°. Curve 
4 is ratio obtained from at the point Z = 80°, 
4=30°. 
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brightnesses and at wavelength 

X0.37 M, computed by the three foregoing meth- 
ods. In the same figure there is the Rayleigh 
extinction coefficient 2* computed for atmosphere 
model CIRA 1965. The dotted line indicates 
the extinction coefficient found for the heights 
60 to 120 km by G. Fiocco and G. Colombo 
(1964) by means of laser observations. The cross 
indicates the extinction coefficient found by A. 
E. Mikirov (1963) by means of rocket measure- 
ments. 

To present the importance of the aerosols 
in the atmospheric scattering, we have computed 
the ratio S/Sj? for different heights from 30 to 


125 km. As is seen from figure 5, the values of 
this ratio are near unity in the height interval 
30 to 80 km and rapidly increase above 80 km. 
It seems that this phenomenon is the result of 
the light scattered by dust which penetrates 
into the Earth's atmosphere from the inter- 
planetary space or arises from meteor ablation. 

The presented data and the results of com- 
paring the twilight data with data obtained by 
other methods permit the conclusion that the 
twilight method offers good possibilities for 
investigation of the optical properties of the 
upper atmosphere and can be used for monitoring 
the dust concentration in the upper atmosphere. 
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